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Aims. We present new spectropolarimetric observations of the chromospheric Hei 10830 A multiplet observed in a filament during 
its phase of activity. 

Methods. The data were recorded with the new Tenerife Infrared Polarimeter (TIP-II) at the German Vacuum Tower Telescope (VTT) 
on 2005 May 18. We inverted the He Stokes profiles using multiple atmospheric components. 

Results. The observed He Stokes profiles display a remarkably wide variety of shapes. Most of the profiles show very broad Stokes 
/ absorptions and complex and spatially variable Stokes V signatures. The inversion of the profiles shows evidence of different 
atmospheric blue- and redshifted components of the He i lines within the resolution element (~ 1 arcsec), with supersonic velocities 
of up to ~ 100 km/s. Up to five different atmospheric components are found in the same profile. We show that even these complex 
profiles can be reliably inverted. 
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1. Introduction 

Filaments are typically elongated dark structures on the solar 
disk, characterized by additional absorption in strong chromo- 
spheric lines, like the Ha line. They are the on-disc counter- 
parts of prominences seen above the solar limb and represent 
dense and cool chromospheric gas suspended in the hot and thin 
corona. The strength and structure of the magnetic field plays 
a key role in keeping the prominence material from flowing 
down to normal chromospheric levels. It is therefore of con- 
siderable interest to determine the magnetic structure associated 
with a prominence. A number of observations of the magnetic 
vector in prominences have b een carried out (e.g., Lerov 1989; 
Casini et al. 2003; Lin et a l.1ll998l; iTruiillo Bueno et al. ibool 



IMerenda et al. »2006i; iKuckein et al. 



2009h . All these investiga- 



I tions, except for the last one, have been restricted to quiescent 
prominences. 

, In this paper we introduce spectropolarimetric observations in 
■ the He 1 10830 A triplet of a filament located in an active region 
' that has been activated by a flare. Filaments also appear as dark 
absorption features in the He i 10830 A line. 
Spectropolarimetry in the He: triplet at 10830 A is an impor- 
tant tool to determine the magnetic field ve ctor in the sola r 
chromosphere (Truiillo Bueno et al. 2002; Sol anki et al. II2003I) . 
It is well suited for the measurement of the magnetic vector 
near the base of the solar corona because these lines are of- 
ten naiTow, nearly optically thin, have a reasona ble effective 
Lande factor (see Table [TJ and are easily observed (iRiiedi et al. I 
|1995h . Measurements of the polarization of the Hei 10830 A 
radiation are also indicated as a useful new tool for t he di- 
agnostics of the mag netic field in fila ments (iLin et a . I 



agnostics or ttie mag netic neld in nlar nents (ILin et a 
iTruiillo Bueno et al. I 12002). Recently, iKuckein et al. 
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studied the vector magnetic field of an active region filament 
by analyzing spectropolarimetric data in the He i 10830 A lines 
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with three different methods, one of these being, as in our 
case, Milne-Eddington inversions. They find the highest field 
strengths measured in filaments so far, around 600-700 G, and 
conclude that strong transverse magnetic fields are present in ac- 
tive region filaments. 

The He 1 10830 A multiplet originates between the atomic levels 
2 1 and 2 ^^"2,1,0- It comprises a component at 10829.091 1 A 
with Ju - (hereafter referred to as Trl), and two components 
at 10830.2501 A with 7„ = 1 (Tr2) and at 10830.3397 A with 
7j, = 2 (Tr3) which are blended at solar chromospheric tempera- 
tures. In the subsequent discussion we call this blended absorp- 
tion the Tr2,3 component. 

We concentrate on describing the observations and the inversion 
of these profiles. The He Stokes profiles in the activated filament 
reveal a remarkable level of complexity, requiring up to five in- 
dependent magnetic components to reproduce. In a following pa- 
per the obtained maps of the magnetic and velocity fields in the 
filament will be presented, critically analyzed and interpreted in 
terms of prominence models. 

2. Observations 

On 2005 May 18 we observed the active region NOAA 10763, 
located at 24° W, 14° S on the solar disk, which corresponds to 
a cosine of the heliocentric angle, fi = 0.9. During the obser- 
vations a flare of GOES class C2.0 erupted in the western part 
of the active region. A filament was present in the active region 
close to and partly overlying the flare ribbons, which was acti- 
vated by the flare. 

The analyzed data were recorded wit h the second genera- 
tion Tenerife Infrared Polarimeter TIP-II dCoUados et al. Il2007h 
mounted at the 70 cm aperture Vacuum Tower Telescope (VTT) 
at the Teide observatory in Tenerife. The TIP-II instrument al- 
lows for a spectral window of 11 A around the He 10830 A 
line to be recorded with a wavelength dispersion of 11 mA 
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Fig. 1. Intensity maps at different wavelengths of the scanned portion of active region NOAA 10763 (24° W, 14° S). Left (a): 
Continuum. Right (b): He i line core of the Tr2,3 component. The four symbols overplotted on the right image refer to the line 
profiles shown in Fig. |2] 



Table 1. Spectral lines contained in the spectral window of the 
TIP-II instrument (10825-10836 A). 



Line 


Wavelength (A) 


Effective Lande factor 


Si I 


10827.09 


1.5 


Hei(Trl) 


10829.09 


2.0 


Cai 


10829.27 


1.00 


Hei(Tr2) 


10830.25 


1.75 


He I (Tr3) 


10830.34 


1.25 


H2O (telluric) 


10832.11 




Cat 


10833.38 


1.00 


H2O (telluric) 


10833.90 




Nai 


10834.85 


1.03 



per pixel. The wide wavelength range turned out to be abso- 
lutely necessary to cover the full range of possible He absorp- 
tion signatures observed during the evolution of the filament. 
As we shall show below in Sect. H] some of the He i profiles 
we observed display absorption at redshifts corresponding to as 
much as 100 km/s. Th is is considerably high er than the highest 
redshifts measured bv Aznar Cuadrado et al. I (I2005.) . whose ex- 
haustive analysis of supersonic downflows seen in the He i lines 
was limited by the spectral range available to the much smaller 
TIP-I detector (Martinez Pillet et al. 1999). 
The spectral window of TIP-II (10825-10836 A) contains pho- 
tospheric Hnes of Sii at 10827.09 A, Cai at 10829.27 and 
10833.38 A and Na i at 10834.85 A, the chromospheric He i mul- 
tiplet and two telluric blends at 10832.11 A and 10833.90 A (see 
Table©. 

The filament was scanned in steps of 0.35" perpendicular to the 
slit orientation (180.51° with respect to the solar N-S direction), 
from 14:38:29 to 15:02:26 UT, providing a map of the region of 
36.5 X 25 Mm2. 

Figure[T]shows Stokes / maps of the observed region obtained in 
the continuum at 10825.8 A (a) and in the core of the Tr2,3 com- 
ponent of the He i triplet (b). The slit is oriented along the y-axis 
and the direction of the scan is along the x-axis. Therefore, along 
the X-axis both spatial and temporal information are mixed. In 
addition to a few small pores. Fig. [T^ clearly shows the granu- 
lation across the whole image, suggesting that the seeing condi- 
tions during the scan were particularly good and stable, allowing 
on average a resolution of 1". The resolution was determined by 
performing a Fast Fourier Transform of the average continuum 
image. Note that the diffraction-limited performance of the VTT 
at this wavelength corresponds to a resolution of 0.39". At some 
locations the He i profile is so broad that no pure continuum is 
present on the detector. These locations appear as dark cloud-like 
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Fig. 2. Observed Stokes / profiles at different locations, indicated 
by the different symbols in Fig.[TJ). Solid line (star symbol): quiet 
Sun profile. Dashed line (cross): Hei emission profile. Dash- 
dotted line (triangle): very broad He i profile indicating multiple 
redshifted atmospheric components. Dotted line (square): deep 
absorption in the He lines, that is in part strongly blueshifted. 
The arrows mark the wavelengths at which images are shown in 

Fig. a 

features in Fig.[T^. In Fig.[TJ) a strong absorption from the fila- 
ment material is clearly visible in the He i lines at 9" < x < 33", 
scanned from ~ 14:43 to 14:53 UT. 

In the last part of the map (33" < x < 49"), at y w 19"-28", 
during the scan of a flare ribbon, the absorption in the He i lines 
is weaker than in quiet Sun regions and the spectra also show 
in part emission profiles (visible as brightenings in the image; 
cf. Lagg et al. 2007). In this region the Ha slit-jaw images that 
were simultaneously recorded at the VTT display flare ribbons 
as well. 

3. Analysis of the Stokes profiles 

In Fig.|2]four observed Stokes / profiles normalized to the con- 
tinuum intensity are plotted as an example of the unusually wide 
range of profiles observed. These spectra and all the spectra we 
show here are binned over four spectral pixels (resulting in a 
wavelength dispersion of 44 mA) and eight spatial pixels (re- 
sulting in a spatial pixel size of 0.7 x 0.7 arcsec^, a single pixel 
being 0.35 x0.18 arcsec-) in order to increase the signal-to-noise 
ratio to ~ 4000. The solid line represents a profile observed in a 
comparatively quiet region of the map (at the location of the as- 
terisk in Fig. [TJi), the dashed spectrum shows a Hei profile in 
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Fig. 3. Example of an observed averaged Stokes vector (dash-dotted profile of Fig. |2]l. The unusually broad Stokes / profile and the 
corresponding signatures of the V profile suggest that the He triplet samples multiple atmospheric components. The vertical dotted 
lines mark the rest wavelengths of the different lines that are present in the observed wavelength range and are identified at the top 
of the figure and listed in Table[T] 



emission observed in a flare ribbon (cross in Fig.[TJ)). Note that 
the Si I line is weakened at this location. The dash-dotted line 
(triangle in Fig. shows a very broad absorption in the He i 
lines observed in the filament. Significant absorption covers the 
wavelength range from 10829 A to 10835 A, with weak absorp- 
tion to the very edge of the detector at 10836 A. The Stokes / 
profile exhibits multiple minima that are not corresponding to 
the rest wavelengths of lines normally found in the Sun's spec- 
trum, nor to telluric absorptions. The multiple absorption dips 
in the profile suggests multiple atmospheric components with a 
shift in the Hei absorption caused by different velocities. For 
the dash-dotted profile these components are nearly at rest or are 
redshifted (see also Fig. [3]). The maximum redshift corresponds 
to a downflow velocity along the line-of-sight of ~ 100 km/s. 
The dotted line (square in Fig. [TJ)) is an example of a particu- 
larly deep absorption in the He i line, which becomes as deep 
as the strong Si i line in the quiet Sun. In addition to this strong 
unshifted absorption there is a broad blueshifted absorption that 
blends with the Si i 10827 A line and extends even bluewards of 
it, revealing rapidly upward propagating chromospheric gas. The 
maximum blueshift corresponds to a velocity along the line-of- 
sight of ~ 60 km/s. As we already pointed out, the shown spectra 
are always binned, but multiple components are also seen in the 
single profiles forming the composite. 

In Fig. [3] we plot the (averaged) full Stokes vector for the dash- 
dotted profile of Fig. |2] (position x = 32", y - 22" in Fig. [TJ. 
The polarized Stokes parameter displaying the highest signal- 
to-noise ratio is Stokes V. It displays a significant signal over 



nearly the entire wavelength range over which Stokes / exhibits 
a strong absorption. The complexity of the Stokes V profile con- 
firms that the broad He i absorption is most likely caused by mul- 
tiple atmospheric components harboring gas that flows at differ- 
ent speeds. It also indicates that at least a part of the downflow- 
ing gas is located in magnetized regions of the solar atmosphere. 
Stokes Q and, to a lesser extent. Stokes U also display a signal 
over a broad wavelength range, although the signal is weak and 
consequently the signal-to-noise ratio is fairly low. In Fig.|3]the 
wavelength positions at rest of the different lines present in the 
observed wavelength range are indicated by the vertical dotted 
lines (see also Table [TJ. 

In order to investigate how the multiple He components are re- 
lated with the motion of the filament material, we show in Fig.|4] 
two intensity maps of the scanned active region at wavelength 
positions -1.08 A (A) and +1.08 A (B) with respect to the rest 
wavelength of the Tr2,3 component of the He lines. The wave- 
lengths at which the images are shown in Fig. |4] correspond to 
up- and downflows of 30 km/s, respectively, and are marked by 
arrows in Fig. |2] From these maps we can see that strong blue- 
or redshifted line profile components are present over practically 
the whole filament, suggesting strong dynamics throughout this 
active filament. These dynamics can be very complex, because 
we note a number of locations at which profiles with strong blue- 
and redshifts coexist within a resolution element. Consequently, 
upflows and downflows are present at the same spatial posi- 
tion, although with this simple analysis alone we cannot rule out 
that these profiles are simply strongly turbulently broadened, or 
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a strong Trl mimics a blues hift. Indeed, a similar anomal ous 
broadening was observed by I Lopez Ariste & Casini I (l2006l) in 
the He 1D3 line in spicules. They found that a distribution of ve- 
locities with FWHM=50 km/s was required to reproduce the ob- 
served broadening. 

In Fig.|5]we present some examples of observed averaged Stokes 
/ (A) and Stokes V (B) profiles at different positions marked by 
the letters a-f in Fig.|4] They stretch approximately from one end 
of the filament to the other and display a gradual change from 
mainly redshifted (a) to mainly blueshifted (f) profiles. All these 
profiles are characterized by multiple atmospheric components 
of the He I lines. In some of them, most prominently in c and 
e, blueshifted and redshifted components coexist. Note also the 
changing ratios of the Si i and the He i Stokes V profile ampli- 
tudes. At some locations, notably b and c, the photospheric and 
chromospheric Stokes V profiles have opposite polarities. This 
indicates that the magnetic structure changes quite dramatically 
between the photosph ere and the chromosphere, which i s not 
unusual for filaments (lLerovl[T989l : iLow & Zhang 1 120021) . We 
also have evidence of different He components showing emis- 
sion and absorption within the same profile (not shown in this 
paper), which are generally co-located with the bright H a flare 
kernels, but do not necessarily appear bright in Fig.[TJ). 

3.1. Stokes inversion 

We analyze the spectropolarimetric observations with the nu- 
merical code HELD£] for the synthesis and inversion of Stokes 
profiles in a Milne-Eddington atmosphere (Lags et al. 2004) ex- 
tended to include the Pasc hen-Back effect jSasso et al. I l2006t 
ISocas-Navarro et al. 1 12004'). HELIX allows for inversions of a 
single spectrum using multiple atmospheric components, so that 
atmospheric parameters like the magnetic field vector and the 
line-of-sight (LOS) velocity are retrieved for each component. 
The free parameters tuned by the code to fit the observed Stokes 
profiles are the magnetic field strength B and direction (inclina- 
tion angle y and azimuthal angle x)^ the line-of-sight velocity 
vlos , the Doppler width AAo, the damping constant F, the ratio 
of the line center to the continuum opacity 770 and the gradient of 
the source fun ction S 1 . For more information on these parame- 
ters we refer to'Lagg et al. ' 1*2004). A range for every parameter 
to be fitted is specified to guarantee that the result of the inver- 
sion stays within the regime of physically useful solutions. 
The inversion includes a scattering polarization correction, a 
treatment of the Hanle effect for a special case. It is only valid at 
disk center and for horizontal fields and it only improves the az- 
imuth retrieval (Truiillo Bueno et al. 2002). Indeed, in the par- 
ticular case of a LOS-direction perpendicular to the solar sur- 
face (disk center) and a magnetic field direction parallel to the 
solar surface (horizontal field) the polarization signal owing to 
the Hanle effect is oriented along the magnetic field direction: 
tan(2;^) = (U/Q). 

The details of the input atmospheres depend on the number of He 
components necessary to reproduce the observed profiles. One 
parameter, the filling factor a,, defines the contribution of a given 
atmospheric component / to the total observed He profile. The 
sum of the filling factors of all atmospheric components within 
a resolution element is required to be unity, 2; ff/ = 1 • The max- 
imum number of He atmospheric components we need to repro- 
duce the profiles in the scan is five (see Fig.|6l). This is consider- 
ably more than in previous studies involving the He line, and we 
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need to demonstrate that these fits are reliable. Additionally, we 
need to determine the maximum amount of information that can 
be reliably gathered from these complex inversions. 
We have tested that we do not employ more components than 
necessary by first fitting each profile with a smaller number of at- 
mospheric components and consequently free parameters. Only 
a significant improvement of the fitness function / (defined in 
Lagg et al. 2004), which gives the goodness of the fit, warrants 
the use of more components or free parameters. Usually, an in- 
sufficient number of free parameters results in a poor fit for parts 
of the observed line profile. 

For some extreme profiles we face the problem of a huge num- 
ber of free parameters to be fitted. The spectral region covered 
by the He absorption signature of strongly broadened or shifted 
profiles also contains four photospheric lines of Si i, Ca i (two 
lines) and Nai and two telluric blends (see Table [TJ. We have to 
fit these lines as well to obtain unbiased fits of the He i absorp- 
tion, which often blends with some or all of these lines. Because 
the photospheric magnetic and velocity vector in general signifi- 
cantly differs from the chromospheric parameters, this increases 
the total number of atmospheric components and requires addi- 
tional line-specific parameters of each line to be determined. For 
simplicity, each spectral line was assigned its own atmospheric 
component. The number of free parameters was limited by cou- 
pling them appropriately (see Sect. 13.2b . To obtain a good fit of 
the strong photospheric Sii line at 10827.09 A, we had to con- 
sider two atmospheric components: a magnetic component and 
a field-free stray-light component. Only this combination could 
satisfactorily reproduce both the line core and wings of the Si i 
line. The two telluric blends result in eight additional free pa- 
rameters, originating from fitting two Voigt profiles to the tel- 
luric lines. The parameters needed to fit the telluric blends could 
be kept in a very narrow range over the whole map, introducing 
no decrease in stability of the convergence of the inversion. 
There are profiles in the scan where the He lines show emis- 
sion features. These profiles were excluded from further anal- 
ysis. Even though the HELIX code can fit emission profiles by 
setting the gradient of the source function 5 1 to a negative num- 
ber, the analysis technique is not optimized to produce reliable 
results. Moreover, in some spatial pixels the He emission pro- 
files coexist with blue- or redshifted He absorptions, making the 
retrieval even more complicated. In the map shown in Fig.|6] the 
pixels containing profiles that were not inverted are black. 



3.2. Parameter coupling 

In order to decrease the number of free parameters and thus to 
enhance the stability of the fitting procedure, we coupled some 
atmospheric parameters between the atmospheric components, 
so that a given atmospheric parameter was forced to have the 
same (but not prescribed) value in all atmospheric components. 
Coupling between parameters of different components was also 
used as a means to impose that all lines from the same atmo- 
spheric layer (in particular the photosphere) must be formed in 
the same atmosphere. We also narrowed the ranges of possible 
solutions for some other parameters, which also helped to ensure 
the uniqueness of the solution. 

In particular, we coupled the magnetic field vector (B, y and x) 
and the gradient of the source function S \ between the four pho- 
tospheric lines, assuming that their values stay approximatively 
constant. This is consistent with the assumption that the four 
lines are formed at similar heights and can be represented by 
a single atmosphere. 
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Fig. 4. Intensity maps of the scanned active region obtained at different wavelength positions indicated by arrows in Fig.|2] Left (A): 
At -30 km/s with respect to the central rest position of the Tr2,3 component of the He triplet. Right (B): At +30 km/s with respect 
to the central rest position of the Tr2,3 component of the He triplet. The different letters refer to the line profiles shown in Fig.|5] 
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Fig. 5. Observed averaged Stokes / (A) and V (B) profiles at the locations marked by letters in Fig.|4] The dashed, vertical lines 
mark the central position of the spectral lines in the wavelength field-of-view of the TIP-II instrument. The (x,y) coordinates in 
arcseconds are written next to each profile. 



Among the photospheric lines, the Si i line is the strongest one 
with the best signal-to-noise ratio and shows clear signatures in 
the Stokes profiles. Changes in the atmospheric parameters will 
have a larger influence on the Stokes profiles of the Si line than 
on the weaker photospheric lines. Therefore, the photospheric 
parameters are mainly determined by the shape and strength of 
the Si Stokes profiles. Then, we made an inversion run for the 
whole map, fitting only the photospheric and the telluric lines. 
From the results of the fit, assuming that these lines are not dis- 
turbed by the filament activity happening at higher atmospheric 
layers, we could select a narrower input range for the Doppler 
width AAo of the photospheric lines and the width and the damp- 
ing of the Voigt profile for the telluric blends. The narrower in- 
put ranges are then used for the inversion done including the He 
lines. This avoids that the code broadens and shifts these lines 
by large amounts to mimic one of the He components. 



The different He components certainly have different LOS veloc- 
ities, but it is not a priori clear if they sample different magnetic 
vectors (and if they do, whether there is sufficient information in 
the measured Stokes parameters to distinguish between different 
magnetic vectors affecting the various Hei absorption compo- 
nents). Therefore, we made several test runs to decide whether 
to couple the magnetic field vector or not (see Sect. HJl. 
We also coupled the gradient of the source function S i between 
the He atmospheric components, and we could also assign a rel- 
atively narrow input range for vlqs of each He component from 
the visual analysis of the observed profiles. Indeed, when the 
velocity separation between the He profiles is distinct enough 
(~ 10 km/s) we can clearly distinguish the positions of the dif- 
ferent minima in the / profile. 
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Fig. 6. Number of He components needed to fit the observed pro- 
files for each pixel of the map. In the map, the pixels containing 
profiles that were not inverted are black (mainly those in which 
the He profile is in emission). 

Table 2. Values of the fitness functions / returned by the fits to 
25 observed averaged profiles. 



Pixel x/y 


78-79 


80-81 


82-83 


84-85 


86-87 


(arcsec) 


(27.3) 


(28.0) 


(28.7) 


(29.4) 


(30.1) 


104-107 


/, = 1.44 


1.42 


1.53 


1.54 


1.29 


(17.7) 


/2 = 1.48 


1.71 


1.75 


1.58 


1.17 




.n = 115 


1.46 


1.29 


1.39 


1.16 


108-111 


/, = 1.95 


1.81 


1.74 


1.86 


1.18 


(18.4) 


/2 = 1.91 


1.55 


1.76 


1.77 


1.04 




/3 = 1.67 


1.48 


1.48 


1.69 


1.08 


112-115 


/, = 1.89 


1.68 


1.73 


1.25 


1.50 


(19.0) 


fi = 2.00 


1.63 


1.60 


1.28 


0.94 




/3 = 1.V8 


1.53 


1.34 


1.23 


1.35 


116-119 


/i = 1.59 


1.76 


1.41 


1.58 


1.56 


(19.7) 


fi = 1.47 


1.97 


1.69 


1.73 


1.52 




/3 = 1.32 


1.80 


1.36 


1.46 


1.40 


120-123 


/, = 1.68 


1.98 


1.51 


1.54 


1.76 


(20.4) 


fi = 1.66 


1.88 


1.51 


1.38 


1.60 




/3 = 1.61 


1.92 


1.53 


1.43 


1.58 



4. Results 

Test runs showed that the best fits to the observed Stokes pro- 
files are obtained by coupling only the magnetic field strength 
B between the He components, leaving the two angles, y and x, 
free. In Fig. |7] we show as an example the results of inversions 
of the observed averaged Stokes profiles (solid lines) at the loca- 
tion (29", 19") of the map in Fig.[T] In Figs.|7ja-d) two fits to the 
observed profiles (dashed and dash-dotted lines) are displayed. 
Both are obtained by allowing five different atmospheric compo- 
nents to describe the He i lines. The dashed lines represent a fit 
obtained by allowing the code to find a different magnetic field 
vector for each atmospheric component used to describe the He 
profiles. For the fit represented by the dash-dotted lines though, 
we coupled the magnetic field strength between the components, 
while leaving the angles describing the magnetic field direction 
uncoupled. The average values of / obtained from 100 inversion 
runs of this particular observed profile are 1.62 for the dashed 
line fit and 1.53 for the dash-dotted line fit. Both fits are good and 
the fitness values are comparable, but for the dash-dotted profile, 
the fit is obtained with 4 free parameters less. The slightly higher 
value of the fitness of the dashed profile is not enough to justify 
the use of so many more free parameters. In Figs.|7je-h) the same 
observed profile (solid line) is plotted as in (a-d). The dotted line 



is obtained by coupling the whole magnetic field vector between 
the five components. The observed Stokes Q, U and V profiles 
cannot be reproduced by this inversion (/ - 0.89). By coupling 
B andx (long dashed lines) we obtain a good fit only for / and V 
i f - 1 .26) and by coupling B and y (dash-dot-dotted lines) we 
can fit /, Q, and U reasonably well, but not V (f - 1.02). In no 
case can we reproduce the Stokes vector as a whole. 
Stokes Q of this observed profile shows a predominantly single- 
lobe signature that suggests resonance scattering polarization. 
The simple scattering polarization correction we used and de- 
scribed in Sect. 13. H is not sufficient to fit the strong Q peak. An 
improved treatment of the scattering polarization may help. 
We carried out the same analysis for another 25 Stokes profiles 
with 3, 4, and 5 He components in a region of 3.5 x 1.8 Mm^ 
around the profile displayed in Fig. |7] The values of the fitness 
function / obtained from three different types of fits (leaving the 
magnetic field free, /i, coupling B, fi, and coupling B 
f^) to the observed profiles are grouped in Table |2] The profiles 
are located at the spatial pixels listed at the top (x) and on the 
left (y) of Table|2] These observed profiles can be reproduced al- 
most equally well either by coupling the magnetic field strength 
between the He components or by leaving the full magnetic field 
vector free. Significantly poorer fits to the full Stokes vectors are 
obtained by coupling B and x- Indeed, < f\ in 88% of the 
cases and /3 < fi in 80% of the cases. f\ is higher than /2 in 
56% of the cases, while in over 40% of the cases /2 > f\ - The 
difference between the two is insignificant and the marginally 
higher value of fx is not enough to justify using a higher number 
of free parameters. 

Another test is based on the analysis of the errors in the retrieval 
of the atmospheric parameters. The convergence of the genetic 
algorithm used for the minimization of the merit function, Pikaia 
(ICharbonneau |[T995h . does not depend on the choice of initial 
guesses for the atmospheric parameters. The convergence to the 
global minimum occurs on a random path after an infinite num- 
ber of iterations. We set an upper limit to the number of iterations 
(700) in a way that the necessary computing time and the stabil- 
ity of the resulting fitness, /, are in good balance. The variations 
of the parameter values resulting from numerous inversion runs 
with an upper limit to the number of iterations can be used as a 
direct measure of the error in the retrieval of this parameter. This 
error contains the intrinsic uncertainty in the parameter retrieval, 
caused by the similarity of Stokes profiles for different parameter 
values, and the error introduced by the inversion process itself. 
The mean retrieved values of the atmospheric parameters with 
their errors, as deduced from 100 inversion runs, are listed in 
Table [3]for each of the profiles shown in Fig.|7] The five values 
per parameter correspond to the different He atmospheric com- 
ponents. As pointed out earlier both the dashed (no coupling) 
and the dash-dotted profiles (only B coupled) provide almost 
equally good fits, but the latter is obtained by an inversion with a 
lower number of free parameters. The errors when B is coupled 
are smaller than without coupling, in particular the errors on B 
and y. This means that with a lower number of free parameters 
the code is more stable. The vios is always well retrieved with a 
significantly lower error than the specified range for the vios of 
any given component (~ 10 km/s). Both the magnetic field an- 
gles, y wdx, can be determined for the individual components. 
The error in the inclination angle, y, is sufficiently large that not 
all components can be clearly distinguished from each other (e. 
g. components 2 and 3 in Table |3]l, but the values obtained are 
in general good enough to deduce the polarity of the magnetic 
field. Regarding the azimuth, the signal-to-noise ratio in the Q 
and U profiles is barely high enough to reliably retrieve infor- 
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mation on it. 

The final three columns of Table [3] list the retrieved values with 
uncertainties for an increasing number of parameters that are 
forced to have the same value in all atmospheric components. 
Irrespective of whether we keep B and y (dash-dot-dotted), B and 
X (long dashes) or B, y and x (dotted) coupled, the fitness / is 
significantly lower than if all parameters are uncoupled, or only 
B is coupled. Furthermore, the spread (uncertainty) in the re- 
trieved parameters tends to be larger in the last three columns, so 
that coupling additional parameters does not lead to an increase 
in the reliability of the retrieved solution, but exactly the oppo- 
site. Furthermore, the retrieved values of the coupled parameters 
are close to the a-weighted averages of the uncoupled values of 
the same parameter. The exception is B, for which the coupled 
value is always lower. Based on the above test (and similar tests 
applied to the Stokes profiles in other pixels) we inverted all the 
observed profiles in the observed map by coupling the magnetic 
field strength between the He components, and leaving x and y 
free (fitness /a in Tabled. 

Examples of the inversions of increasingly complex Hei 
10830 A Stokes /, Q, U, and V profiles ai-e shown in Figs. [Mil 
In each plot, the black lines are the averaged observed profiles, 
while the red lines are the best fits to the observations obtained 
considering one, two, three, four, and five atmospheric compo- 
nents of the He I lines. These are the minimum numbers of at- 
mospheric components that are needed to obtain a satisfactory 
fit, both as judged by eye and from a decrease of / found by 
adding a further component. Also plotted are the different He 
components (magenta, light blue, yellow, green, and dark gray 
lines), whose sum gives the red lines. The plotted contribution 
of each component is multiplied by its filling factor. In Table |4] 
we give the mean retrieved values of the atmospheric parameters 
B, y, X, vlos and a, together with their errors. The values listed 
in Table|4]are averages of the parameter values returned from 50 
inversion runs, and the uncertainties reflect the spread (standard 
deviation) of the returned values. When multiple-component fits 
to the He line are required, such as those displayed in Figs.l9HT2l 
the parameters retrieved for each component are given. Note that 
the photospheric lines are simultaneously fitted by two atmo- 
spheric components (one magnetic, one field free), as described 
in Sect. 13.11 For each fit we also give the value of the fitness 
function in the Table, / + A/. 

Figure|8]displays a one-He-component fit to the observed Stokes 
profiles. Similar quiet profiles are found in the region of the scan 
outside the portion with the filament or other activity (flare). In 
Table |4] (first table) the mean values and the errors of the re- 
trieved atmospheric parameters after 50 inversion runs are re- 
ported. Even with this simple profile we have difficulties in re- 
trieving the values of the magnetic field angles, y and x, in par- 
ticular X with high accuracy. This is mainly because of the low 
signal-to-noise ratio of the observed Q and U profiles and to a 
smaller extent because of the already high number of free param- 
eters to be fitted (the inversion includes the He, the photospheric 
and the telluric lines). Figures|9]and[T0]show a two- and a three- 
component fit to other, more complex observed profiles. The two 
profiles, observed close in position and time, are somewhat sim- 
ilar, with the difference of a deeper He blueshifted absorption 
in the latter profile. This additional absorption requires an ad- 
ditional He component to be fitted. The presence of another He 
component in the profile of Fig. [TO]is also evident from the dif- 
ferent shape of the negative lobe in the Si i V profile. In both 
profiles we find strongly supersonic blueshifted He components 
(speeds in excess of -50 km/s) coexisting with a He component 
nearby at rest. The fit to the observed profile shown in Fig. [TT] 



requires an additional fourth redshifted He component. The fi- 
nal fit shown in Fig. fT2] requires five atmospheric components. 
Only few other profiles (see Fig. |6]l need five components to be 
reproduced, and these are all located around one end point of the 
observed filament. This is also the location where we measure 
the highest downflow velocities (~ 100 km/s). It is particularly 
gratifying to see that both the fitness / and the uncertainties of 
the retrieved parameters appear to be nearly independent of the 
number of atmospheric components. This is probably because 
profiles with more components tend to have a larger equivalent 
width and individual components are well separated in wave- 
length. 



5. Discussion and conclusions 

We focus on the analysis and the inversion of Stokes profiles 
of the Hei 10830 A triplet observed in an active region fil- 
ament during its phase of activity. We showed that even the 
most complex observed profiles can be well reproduced by 
multi-component atmospheres. Although we cannot completely 
exclude that the profiles could be simply strongly turbulently 
broadened, the fits based on multiple components are of clearly 
better quality (see below). We tested the response of the numeri- 
cal code HELIX for the inversion of Stokes profiles to work with 
a high number of atmospheric components and free parameters. 
This was necessary to fit not only the different observed He com- 
ponents, but also the photospheric and telluric lines in the wave- 
length range. We were able to retrieve the atmospheric parame- 
ters for the individual He components at different levels of ac- 
curacy. The values of the magnetic field strength, B, retrieved 
from the fit presented here are in the range 100 - 25 G. Th ese 
values are lower than those found by i Kuckein et al. I (l2009h for 
another active prominence. This difference coul d be caused by 
the different methods employed, because unlike iKuckein et al. I 
(2009) we inverted our profiles including a scattering polariza- 
tion correction (see Sec. 13. Il l and the magnetic field strength we 
retrieved is a result of a coupling between many He components. 
The difference could also be merely a reflection of the different 
strengths of the two filaments' magnetic fields, or the analyzed 
profiles possibly do not sample the highest field strength region 
in the observed filament. The v^os is well retrieved for all pro- 
files with a small error. The inclination angle y is retrieved with a 
somewhat larger error, but the inclinations of the different com- 
ponents of the magnetic field can often be distinguished. The 
error bars on the azimuth angle, x^ are instead quite large, and 
it is often impossible to distinguish between the azimuths of the 
different magnetic components. This is mainly because of the 
complexity and the noise in the Q and U Stokes profiles. The 
information in the observations for retrieving the azimuth is lim- 
ited. We therefore refrain from conclusions about the azimuthal 
direction of the magnetic field in the filament. There is also in- 
sufficient information in the profiles to distinguish between the 
field strength B of the different components. Note that because 
Stokes V is generally much stronger than Q and U, B cos y is 
probably the most reliably determined quantity, although we do 
not explicitly tabulate it. 

The analysis of the observed polarization of the Hei 10830 A 
multiplet in the filament, carried out by inverting the Stokes pro- 
files, reveals different unresolved atmospheric components of 
the He lines, coexisting within each spatial resolution element 
(~ 1 arcsec). The components are distinguished by their Doppler 
shifts, which generally differ by more than the sound speed. For 
more complex analyzed profiles we also tried to find a solution 
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Fig. 8. Observed spatially averaged and inverted Stokes /, Q, U and V profiles at the location (5", 13") of the map in Fig.[T] The 
solid black lines are the observed profiles, while the red lines are the fits to all observed lines obtained imposing one atmospheric 
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Parameter/Fit-profile 


Dashed Dash-dotted Dash-dot-dotted Long dashes Dotted 
No coupling B coupled B,y coupled B,x coupled B,y,x coupled 


B±AB (G) 


478 ± 330 228 ± 53 250 + 134 234 ± 55 222 + 123 
264 ± 215 
403 ± 285 
307 ± 159 
632 ± 376 


7 ± Ay (°) 


47 ±25 33 ±22 122 ± 23 31 ±26 124 ± 22 
122 ±22 110 ±4 110 ±5 
108 ± 24 102 ±7 102 ± 7 
131 ±21 141 ±17 143 ±17 
141 ±22 170 ±10 169 ±12 


X±^X (°) 


-29 ±34 -18 ±48 -30 ±43 30 ± 42 27 ± 44 
46 ±70 52 ± 62 60 ± 53 
58 ±28 48 ± 32 50 ± 32 

7 ±16 0±13 1±13 

8 ± 33 -2 ± 37 -3 ± 37 


Vlos ± Avlos (km/s) 


-38 ±3 -39 ±2 -31 ±2 -38 ± 3 -31 ±2 
1±1 1.1 ±0.7 0.9 ±1.0 1±1 0.9 ±1.0 
21 ±1 21 ±1 21 ±1 21 ± 1 21 ± 1 
53 ± 1 53 ± 1 53 ± 1 53 ± 1 53 ± 1 
68 ± 2 69 ± 2 68 ±2 69 ± 2 68 ± 2 


a ± Aa 


0.10 ±0.03 0.10 ± 0.02 0.09 ±0.02 0.10 ±0.03 0.10 ±0.03 
0.33 ±0.06 0.32 ±0.04 0.29 ±0.04 0.31 ±0.04 0.27 ± 0.05 
0.20 ±0.05 0.19 ±0.03 0.19 ±0.04 0.21 ±0.03 0.21 ±0.03 
0.25 ±0.06 0.24 ±0.03 0.27 ±0.04 0.23 ± 0.04 0.26 ± 0.04 
0.12 ±0.05 0.14 ±0.03 0.16 ±0.04 0.15 ±0.03 0.16 ±0.03 




1.62 ±0.19 1.53 ±0.16 1.02 ±0.06 1.26 ±0.09 0.89 ±0.04 



considering a broadened profile with some emission contribu- 
tions distributed at the appropriate places, but the fits to the ob- 
served profiles were not as good as the ones we present in the 
paper The main problem is reproducing the / and V profiles at 
the same time. 

In this active filament we find profiles requiring up to five 
atmospheric components for a reasonable fit. Multiple un- 
resolved magnetic components are found not only in fila- 
ments, but also above pores and elsewhere in active regions 
(lAznar Cuadrado et al. I l2005t iLagg et al. I l2007h . but so far in 
ot her data sets never more than three components were needed. 
As'Lagg et alj (l2007h pointed out, the geometrical interpretation 
of these multiple downflow components is not clear-cut in the 
case of a line such as He i 10830 A, which is generally optically 
thin. Thus the different components could be located at differ- 
ent horizontal positions within the resolution element, implying 
considerable fine structure in the filament. This fine structure in 
fil aments at a sub-arc s scale has been reported e arlier See, e. 
g., iTandberg-Hanssen 1 ([T99l or lHeinzelldlOOTi) . for reviews. 
These fine structures show a random motion with velocities of 
5-10 km/s, as deduced mainly from observations in the H a line. 
An alternative explanation is that the various components are lo- 
cated at different heights along the LOS. It may be argued that 
the He: line gets very strong in the filament and can hardly be 
considered to be optically thin, so that one should not be able to 
see through different layers so easily. We note, however, that the 
line becomes optically thick only at wavelengths at which the 
absorption is large, whereas it remains optically thin at neigh- 
boring wavelengths. Consequently, if the different layers absorb 
at wavelengths that are shifted by more than a Doppler width 
relative to each other, this explanation remains possible even if 
individual layers produce optically thick absorption. 
Because in many cases the profiles have some optical thickness 
(they are formed over a range of heights), it is possible that gra- 



dients in velocity along the LOS can affect their shape. The in- 
fluence of these gradients cannot be treated in the simple model 
employed here. A more general treatment with gradients might 
lead to good fits with a reduced number of components. 
The inversions support the idea that the He components corre- 
spond to plasma trapped in different magnetic field lines, which 
may well be pointing in different directions along the LOS. 
Indeed, the coupling of the magnetic field vector between the 
individual components results in worse fits to the observed pro- 
files, which supports this interpretation. Note though that some 
He profiles do not show a significant signal in Stokes V. They 
describe intensity features without a polarization counterpart. 
Multiple magnetic components are fairly common in most parts 
of the observed filament and are often associated with strong 
blue- or redshifts corresponding to supersonic velocities. We 
measured downflow velocities of up to 100 km/s and upflows 
of up to 60 km/s along the line-of-sight. These supersonic up- 
and downflows always coexist with a He atmospheric compo- 
nent almost at rest (-10 < vios < +10 km/s) within the same 
resolution element. Sometimes strong blue- and redshifted com- 
ponents are found to coexist in the same profile. Velocities of 
filament fine structures reported in the literature are low com- 
pared with these observed supersonic down- and upflows. 
Strong motions and velocities in active region filaments were 
previously observed in the He I 10830 A lines. Observations 
of an erupting active region filament presented by Penn . (20dol) 
revealed the He I absorption line blue-shifted to velocities be- 
tween 200 and 300 km/s. In the same observation, the filament 
also showed internal motions with multiple Doppler components 
shifted b y +25 k r n/s. R edshifts of 30-60 km/s were reported 
by P enn & Kuhn I (119951) in an active filament during the decay 
phase of a flare in the He I 10830 absorption lines. We found 
the highest redshifts so far measured with the He i 10830 A line. 
In this respect, we emphasize the importance of the TIP-II in- 
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Fig. 9. Observed and inverted Stokes /, Q, U, and V profiles at the location (27", 9") of the map in Fig. [T] The solid black lines are 
the observed profiles, while the red lines are the fits obtained imposing two different atmospheric components to describe the He I 
lines (the profiles of the individual components are given by magenta and light blue lines). 



strument, which was used for the first time during this obser- 
vation campaign. We were able to observe the full range of He 
absorption signatures and measure the highest downflow veloc- 
ities thanks to its wider range of wavelengths compared to the 
original TIP instrument. 
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Fig. 10. Observed and inverted Stokes /, Q, U, and V profiles at the location (29", 13") of the map in Fig. [T] The solid black lines 
are the observed profiles, while the red lines are the fits obtained by allowing three diff'erent atmospheric components to describe 
the He I lines (represented by magenta, light blue, and yellow lines). 
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Fig. 11. Observed and inverted Stokes /, Q, U, and V profiles at the location (17", 16") of the map in Fig. [T] The solid black lines 
are the observed profiles, while the red lines are the fits involving four different atmospheric components to describe the He I lines 
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Fig. 12. Observed and inverted Stokes /, Q, U, and V profiles at the location (32", 21") of the map in Fig. [T] The solid black lines 
are the observed profiles, while the red lines are the fits involving five diff'erent atmospheric components to describe the He I lines 
(magenta, light blue, yellow, green, and dark gray lines). 

Table 4. Mean values of the atmospheric parameters with their errors returned by the fits in Figs.[8l-fT2l 



fitness / = 2.20 ± 0.05 



fitness /= 1.72 ±0.07 



Parameter 


Component 1 


Parameter 


Component 1 


Component 2 


B(G) 


164 + 88 


B(G) 


141 ± 82" 


141 ± 82" 


r(°) 


121 ±20 


yC) 


82 ±5 


40 ±26 




11 ±51 


xC) 


36 ±5 


38 ±40 


Vlos (km/s) 


3.0 ±0.1 


Vlos (km/s) 


-52 ± 1 


1.8 ±0.2 


a 


1.00 


a 


0.35 ± 0.05 


0.65 ± 0.05 





fitness / = 


2.35 ± 0.06 




Parameter 


Component 1 


Component 2 


Component 3 


S(G) 


135 ±45" 


135 ±45" 


135 ±45" 


r(°) 


83 ±5 


124 ± 19 


27± 19 


xC) 


-10±2 


-16± 12 


4±7 


Vlos (km/s) 


-53 ± 1 


-30 ±2 


1.9 ±0.2 


a 


0.32 ± 0.04 


0.21 ± 0.04 


0.47 ± 0.03 




fitness / = 


1.90 ±0.08 





Parameter Component 1 Component 2 Component 3 Component 4 



B(G) 


106 ± 46" 


106 ± 46" 


106 ± 46" 


106 ± 46" 


7(1 


80 ±7 


91 ±5 


103 ±5 


160 ± 18 


xC) 


-17 ±30 


-17 ±30 


1 ± 19 


35± 11 


Vlos (km/s) 


-52 ±2 


-18± 1 


1 ±2 


19 ± 1 


a 


0.17 ±0.04 


0.23 ± 0.05 


0.36 ± 0.04 


0.24 ± 0.04 



fitness /= 1.79 ±0.11 



Parameter Component 1 Component 2 Component 3 Component 4 Component 5 



B(G) 


244 ± 55" 


244 ± 55" 


244 ± 55" 


244 ± 55" 


244 ± 55" 


rO 


11 ± 11 


126 ± 10 


89 ±6 


154 ± 17 


154 ± 18 


xC) 


15 ±47 


15 ±40 


12 ±40 


-19±9 


3 ±65 


Vlos (km/s) 


-25 ± 1 


3.6 ±0.5 


25 ± 1 


67.3 ± 0.6 


100 ± 1 


a 


0.08 ± 0.01 


0.28 ± 0.04 


0.23 ± 0.04 


0.29 ± 0.03 


0.11 ±0.03 



" The B values in all components were forced to be equal (in each table). 



